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The quark-hadron phase transition at finite baryon chemical potential is investigated in 
the extended Nambu-Jona-Lasinio model in which the scalar-vector eight-point interaction 
is included holding the chiral symmetry. By comparing a pressure of the symmetric nu- 
clear matter with that of the quark matter, the realized phase at given baryon density is 
determined. As a result, the first order quark-hadron phase transition is derived where the 
deconfinement transition occurs after the chiral symmetry restoration in this model. 

^! §1. Introduction 

^ ■ One of the recent interests in the systems governed by the quantum chromody- 

namics (QCD) may be to understand the phase diagram and/or the phase transitions 
between various phases in the quark-gluon and the hadronic matters. As for the chi- 
^ ■ ral phase transition, many works have been done in the various effective models of 

I QCD at finite temperature and/or density. However, to obtain definite results of the 

' quark-hadron phase transition at finite baryon density is still difficult because of the 

■ color confinement on the side of the hadronic phase. Of course, the study based on 
the lattice QCD is also difficult at finite density at present. 

' For the symmetric nuclear matter, it is important to describe the property of 

■ nuclear saturation. The Walecka model^' has succeeded in describing the saturation 
property of symmetric nuclear matter as a relativistic system. The success is mainly 
due to the cancellation of the large scalar and vector potential which are derived by 

^ . the a meson and u meson exchanges between nucleons, respectively. In this model, 

^ I the nucleon is treated as not a composit but a fundamental particle. Although this 

model has given many successful results for nuclei and nuclear matter, this model at 
first stage has no chiral symmetry which plays an important role in QCD. 

The celebrated Nambu-Jona-Lasinio (NJL) modeP' gives many important re- 
sults for hadronic worlci^ based on the concenpts of the chiral symmetry and the 
dynamical chiral symmetry breaking. This model has been applied to the investiga- 
tion of the dense quark matter.'^!' Also, by using this model, the stability of nuclear 
matter, as well as quark matter, was investigated in which the nucleon is constructed 
from the viewpoint of quark-diquark picture.'S' On the other hand, it is known that, 
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if the nucleoli field is regarded as a fundamental fermion field, not composite one, 
the nuclear saturation property can not be reproduced starting from the original 
NJL Lagrangian. However, if the scalar-vector and isoscalar-vector eight-point in- 
teractions are introduced holding the chiral symmetry in the original NJL model, 
the nuclear saturation property is well reproduced^- where the nucleon is treated as 
a fundamental fermion. 

In this paper, paying an attention to the chiral symmetry, the NJL-type model is 
adopted in both the nuclear and quark matters. As for the nuclear matter, we adopt 
the extended NJL model with the eight-point interactions and regard the nucleon 
field as a fundamental fermion field with the number of color, Nc, being one-^^® 
Also, we adopt the extended NJL model with Nc = 2> for the quark matter. Thus, the 
quark-hadron phase transition is investigated by using the same-type model with the 
chiral symmetry for nuclear and quark matters in a unified way, the model parameters 
and the number of color of which are different. As the first step to invetigate the 
quark-hadron phase transition, in this paper, we consider the symmetric nuclear 
matter and free quark phase without quark-pair correlation. Namely, we do not take 
into account of the color superconducting phase in the quark phase at finite density 
and at zero temperature.lS' This phase may exist at finite density, but, in this paper, 
the investigation is rest at future. Also, we only consider the symmetric nuclear 
matter, while it is interesting to the study of the neutron matter which leads to the 
understanding of the physics of neutron star.^^^ 

This paper, is organized as follows: In the next section, we introduce the ex- 
tended NJL model at finite temperature and baryon chemical potential for nuclear 
and quark matters. The mean field approximation in this model is given. The derived 
results are identical to those derived by the minimization of the thermodynamical 
potential density. In §3, the method to determine the model parameters is presented 
for nuclear matter and quark matter, respectively. In §4, the numerical results are 
given. In §5, the quark-hadron phase transition is described in this model. The last 
section is devoted to a summary and concluding remarks. 

§2. Extended NJL model for nuclear and quark matter at finite 

temperature and density 

In this section, the same models for the nuclear and quark matters are consid- 
ered in which parameters are different for nuclear and quark matters. One possibe 
model may be a Nambu-Jona-Lasinio (NJL) model^' for nucleons and/or quarks. 
The NJL model is useful because it implements in a simple way the important chiral 
symmetry. As for the nuclear matter, the property of nuclear saturation must be re- 
alized. Then, as one of the possible models, we give the NJL-type model in which the 
fermions interact with each other through the original NJL-type four-point intarac- 
tion plus four point vector-vector and eight-point scalar-vector and isoscalar-vector 
interactions. This model is called an extended NJL model in this paper. 
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2.1. Extended NJL model for nuclear and quark matters at finite density and tem- 
perature 

In this paper, we start with the following Lagrangian densitj©''^^ for nuclear 
matter (i = N) and qurak matter (i = q): 

= Ip.ij^d^^^ + Gl{{^,^l^,f + (v;,i75rV^,)')_- G'Mri^,){^il,.i^i) 

-GUi^.i^if + (Vii75TV'i)')(V'.7^V'*)(V'i7MV'0 • (2-1) 

Here, represents nucleon field {i = N) or quark field {i = q). The first two terms 
give the original NJL Lagrangian density. However, the nuclear matter saturation 
properties are not reproduced if only these two terms are taken into account. Thus, 
we introduce other two terms, following Refl6]). Namely, a vector- vector repulsive 
term, whose interaction strength is represented by G^, and a vector-scalar coupling 
term, whose interaction strength is represented by G*^, are introduced. This model 
is nonrenormalizable. Thus, it is necessary to introduce the cutoff parameter Ai. We 
adopt a three- momentum cutoff scheme. 

We apply the mean field approximation to the above Lagrangian density. By 
mean field approximation we mean that we replace bi-linear quantities in the fermion 
fields, such as ipFi/j, by {tpFxp) + {ipFip — {ipFip)), and keep only linear terms in the 
fluctuation {ipFip — {ipFip)). Here, the symbol (• • • ) denotes the expectation value 
or thermal average. Only the expectation values of ip^i^j^d^ipi and ipiipi survive at 
finite temperature and density. The mean field Lagrangian density Cf^^ and the 
mean field Hamiltonian density Hf^^ are easily obtained as 

^MF ^ . vV'i + mitpitpi + Jliipi-f^ipi - Ci , 

Q = -G\{(^M? + G;#,7V.f + GL#.V'.f#i7V.f , (2-2) 
where we define m-j and Jli as 

m, = -2G^#,V.» + 2G:,#,V'i))(%7V^f , (2-3) 
/I, = 2Gt((V,7V^)> + 2GL((V'iV'if ((V'^7V.)) (2-4) 

for the nuclear matter {i = N) and the quark matter (i = q), respectively. Here, the 
symbol ((•••)) denotes the thermal average which is given latter, while in this paper 
we only treat the zero temperature system. Also, the symbol (• • • ) will be used for 
the expectation values at zero temperature. 

We deal with a finite density system in this paper. In order to calculate the 
physical quantities at finite density, we introduce the chemical potential : 

= -itpil ■ Vipi + niiipiil^i - iiltpi-f^tpi - d , (2-5) 
where the effective chemical potential ^[ is defined as 

Ml = - = - [2G;#,7V*» + 2G^,#,V^»'#^7Vi»] ■ (2-6) 
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Under this Hamiltonian density, we can calculate physical quantities for nuclear 
matter and/or quark matter. 

The expectation values at zero temperature can be compactly expressed as 

(V^^V'^) = -/ T^Tr(^5.(p)) = -ANlN}m, [ f \ ^ \ . , 
J ^(zvrj^ ' J i[ZT:pp'^ — mi — i€ 

.,0„,, \_ I' A rp,/n-o/..NN A^n^ri f P° 



J i[2tt)* •'J lyiTipp'^ — m^—ie 

mi • Pm = -[ vT^Tr(7 • piS.ip)) = -mlNjm, [ .f\^ , . , (2-7) 

where iSi{p) = l/{p^'j^ — rrii — ie) is the fermion propagator in which po is replaced 
into Po + • We introduce Nj and A*"* which represent the numbers of flavor and 
color. For nuclear matter, = 2 and A'^ = 1 are adopted, and for quark matter 
Nj = 2 and Nc = 3 are adopted. Further, we can use imaginary time formalism to 
deal with the system under consideration at finite temperature. As is well known, 
the time component of four- momentum after Wick's rotation, p4, can be replaced 
by Matsubara's frequency a;„ = (2n + l)7rT where T = 1/(3 is temperature and 
the integral J dpQ /{2tt) is also replaced into the Matsubara sumPJi iuj^ and 

/ (ipo/(27r) — ^^X^^-oo- ^ result, we can calculate the above derived quantities 
at zero temperature with the degeneracy factor i^i = 2NjN^ : 

^^V'.)) = / ^ (nV - nl) , (2-8) 



^^l'^^)) = ^i I (^«+"-) ' (2-9) 

. p)^^)) = [ fP K-nL) , (2-10) 



(2-11) 



where n!|_ is the fermion number distribution functions defined as 



Here, the contribution of the occupied negative energy states should be eliminated 
from the nucleon and/or quark number density itself in Eq. (j2-9p . Therefore, we will 
replace n}_ by n}_ - 1 in Eq.dMD- Thus, Eqs.dMD, dMI), ([MD and ([231) with (f2TT]l 
constitute a set of the self-consistent equations. 

2.2. Thermodynamical potential density 

Equation (|2-3|) with (|2-8|) and (|2-9|) gives a self-consistent equation for nii. This 
is nothing but the so-called gap equation. This gap equation ()2-3p and the fermion 
number distribution functions (j2-lip presented in the previous subsection are ob- 
tained by minimizing the thermodynamical potential density coi in the mean field 



First Order Quark-Hadron Phase Transition in the NJL-Type Model 



5 



approximation. Here, uji is defined as 



u;^ = {{nr))-^^^m))-^{{s^)) 



(2-12) 



where 



(2-13) 
(2-14) 



(2^ 



Inn^ + (1 - n!,_) ln(l - n^) 

+n'_ Inn'l + (1 - n'_) ln(l - n'_) 



(2-15) 



for nuclear matter [i = N) and quark matter {i = q), respectively. Here the ther- 
modynamical averages ((O)) = TrC'e"^^^* -/iiM)/Xre~^(^' -/^iM) j^g^yg ^j^g same 
forms as Eqs. (|2-8p ^ ()2-10p . By minimizing with respect to rrij, n!,_ and nL, we 
get the following equations : 



:{n\-nl) 



X [-mi - 2GlU,iP,)) + 2GL#.7>if P.V'^))] = , 
1 



(2-16) 



duji 



p2 _|_ ^2 



|p2 ^ ^^(_2G*^((v^,v^)) + 2GlAi^a'^,fU^i^^)))} 



-(/X, - 2G;((V',7>.)> - 2GlAi^a'i^i)){{i;,^l;,f) + ^ In 



/3 1 - 



, 

(2-17) 



dn'[_ 



p2 _|_ 



/3 1 - n! 



. 

(2-18) 



Thus, we obtain the gap equation (j2-3p from (|2-16p again, and then, the fermion 
number distribution functions n!|_ in (|2-lip are also obtained from Eqs. (|2-17p and 
(|2-18p with the effective chemical potential ()2-6p . 



§3. Determination of model parameters 



In this section, we give the method to determine the model parameters. 
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3.1. For nuclear matter 

In this subsection, we give the prescription to fix the model parameters whose 
numerical values are presented in the next section. In the model for nuclear matter 
given in the previous section, we have four parameters: G^, G^, and the three 
momentum cutoff A]\f. We can fix these parameters so as to reproduce the vacuum 
and saturation properties for nucleon and nuclear matter at zero temperature. At 
zero temperature, the fermion number distribution function given in Eq. ()2-lip is 
reduced into 

= e{^i^■^ - ^p^+m%) , = 1 , (3-1) 

where 6 is the Heaviside step function. Here, the contribution of the anti-nucleon 
for nuclear matter density in Eq. (j2-9p for i = is subtracted as was mentioned in 
the previous section, namely, should be replaced to — 1. Thus, the nucleon 
number density ()2-9p is calculated as 



where we have introduced the Fermi momentum pp = y /i^^ — for nuclear 
matter. Then, the gap equation ()2-3p with ()2-8p can be expressed as 



rriN = -2Gf (l - (V'^Vtv) , 

Wn^^n) = / d\p\ , = . (3-3) 



It is seen from ()3-3p that the scalar-vector coupling term in the Lagrangian density 
(I24|) gives the density dependent coupling G^{pn) (= Gf (1 - G^^/G^ ■ p%)) in 
the gap equation of the original NJL model Lagrangian with G^ = G^, = 0. Also, 
{'4'n{i ' p)'4'n) is calculated as 

&N{l-P)i^N) = d\p\ ^^1— . (3-4) 



Further, the energy density per single nucleon at finite baryon density and zero 
temperature is easily evaluated as 

Wn{pn) = ruNiPN = 0) , (3-5) 

PN 

where 

(nDiPN) = ii^Nh ■ P)M - Gf (1 - GfjGf • p%){ip^M' + G^p% . (3-6) 



Here, we have denoted the density dependence explicitly. 
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From ()3-3p . we fix the nucleon mass at zero and normal nuclear density as 
iTT-NiPN = 0) = 939 MeV and mN{p^) = 0.6mN{pN = 0), respectively. Here, 

= 0.17 fm~^ is a normal nuclear density and mj\f{p^) is given as a typical 
effective mass value at the normal nuclear density. The reason why we take the 
effective nucleon mass mN{p%) = 0.6mN{pN = 0) is that the quark mass in the 
original NJL model with Giv = Gv = can be numerically calculated as 187 MeV 
at p^. Thus, the nucleon mass at p^ is approximately evaluated as 561 = 3 x 187 
MeV which is about O.Gm^ (~ 563 MeV). From (|3-5p . we fix the parameters so as 
to reproduce the nuclear matter saturation properties, namely, Wn should give the 
minimum value —15 MeV at normal nuclear density p^. Then, our model-parameters 
are determined by these four conditions: nucleon mass at zero and normal nuclear 
density and the saturation properties. As will be shown in §§4.1., the nucleon mass at 
normal nuclear density gives an influence to the imcompressibility of nuclear matter 
at normal nuclear density. Therefore, it is possible that the imcompressibility of 
nuclear matter at normal nuclear density is adopted as an input parameter instead 
of the nucleon mass at normal nuclear density. However, in this paper, we fix the 
value of the nucleon mass at the normal nuclear density derived by the simple quark 
model which is implemented by the quark NJL model. 

3.2. For quark matter 

For quark matter, we use same model Lagrangian density (|'2-1|) as that used 
in the nuclear matter with different model parameters. For quark matter, the pa- 
rameters Gv and Giv are usually set equal to 0, because, in nuclear matter, the 
corresponding parameters were introduced so as to reproduce the nuclear matter 
saturation properties. Thus, we put = 0. However, we retain Giv which appears 
in the gap equation (|2-3p at finite density. Then, the scalar-vector and isoscalar- 
vector interaction term with Giv gives an influence to the chiral phase transition 
at finite density. Namely, Giv may be used to tune the chiral phase transition and 
the slope of the equation of state at high densities, while Giv is treated as a free 
parameter in this paper. 

Then, there are three model parameters, namely, Gf , the three momentum cutoff 
Aq and Giv As for the Gf and Ag, these two parameters are determined by giving 
the dynamical quark mass mg, which is obtained by solving the gap equation ()2-3p 
in the vacuum, and the pion decay constant fj^: 



Jo 



2 



+ ml 



fl = ^ rf|pl7-r^-^- (3-7) 



Then, and Aq are taken so as to reproduce ruq = 313 MeV and /tt = 93 MeV in 
the vacuum. 

Here, Glv is taken as a free parameter in this paper. The change of G'iv leads 
to the change of the chiral phase transition point as will be given later. 

The physical quantities at finite density and temperature can be calculated sim- 
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ilar to the case of the nuclear matter. Then, rriAr, /i^ and z^at should read ruq, //^ and 
Vq = 2NjNc with Nj = 2 and Nc = 3, respectively, and Gt = 0. Here, fiq is a quark 
chemical potential. Of course, in the previous expressions, should be replaced to 
n'L — 1 in the expression of the quark number density. 



§4. Numerical Results 



4.1. Nuclear Matter 



For nuclear matter, we take the parameters so as to reproduce the nucleon mass 
at normal nuclear density as mjv(p^) = 0.6mjv(pjv = 0). Other conditions are 
fixed as follows: The nucleon mass in vacuum, m]\j = 939 MeV, and the saturation 
properties of infinite nuclear matter, Wm = — 15 MeV at the normal nuclear density 

= 0.17 fm~^. Then, the parameter set is summarized in Table I. 

In these model parameters, the incompressibility of the nuclear matter at the 
normal nuclear density is numerically evaluated as 



^ = ^Pn T2 



Pn=P% 



260 [MeV] . (4-1) 



Thus, a rather reasonable value is obtained. In Fig{Tl the energy density per single 
nucleon in Eq. ()3-5p is depicted as a function of the nuclear matter density divided 
by the normal nuclear matter density p^, namely, Pn/p%- 

However, the momentum cutoff A^- is rather small. This leads to the question 
of the applicability of this model. The momentum of nucleon should be smaller than 
the value of three-momentum cutoff. We thus do not take into account of the scalar 
excitation of nucleon and antinucleon whose mass is the twice of nucleon mass. 

If we take the nucleon mass value at the normal nuclear density as mj\f{p^) = 
0.75mN{pN = 0), then the parameters are obtained as An = 441.2 MeV, GfA% = 
16.7483, G^^A% = -1958.12 and G^ A% = 16.9872. Then, the imcompressibility 
is obtained as K = 286 MeV. Thus, the equation of state becomes stiff at normal 
nuclear density. 

4.2. Quark Matter 

For quark matter, we set Gv = in Eq. (|2Tp as was already mentioned in the 
preceeding section. As for the model parameters, the coupling constant Gf and the 



Table I. The values of model parameters for nuclear matter. 



An 


377.8 [MeV] 




19.2596 




-1069.89 




17.9824 



Table II. The values of model parameters for quark matter. 





653.961 [MeV] 


GUI 


2.13922 
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binding energy per nucleon [MeV] 




nuclear matter density divided by normal nuclear density 

Fig. 1. Nuclear matter saturation properties with miv(p5v) = 0.6mjv(pjv = 0): Parameters are 
summarized in Table I. Then, the imcompressibility is A' = 260 [MeV]. 



quark mass [MeV] 
mq 




315 320 325 330 p,^ 

quark chemical potential [MeV] 



Fig. 2. The solutions of the gap equations for quark mass are depicted as a function of the quark 
chemical potential fiq with Giy — for the original NJL model. The vertical axis represents the 
dynamical quark mass niq and the horizontal axis represents the quark chemical potential fiq. 



three momentum cutoff Aq are determined from tlie dynamical quark mass and tlie 
pion decay constant in tlie vacuum in Eq. (|3-7p . Tfie parameters are summarized in 
Table II. 

4.2.1. Giv = 

First, we use the original NJL model Lagrangian density without the vector- 
scalar (Giv) interactions, namely, we put Giv = 0. We solve the gap equation in 
Eq. (|2-3p with Giv = under a finite quark chemical potential at zero temperature. 
Then, the dynamical quark mass niq is obtained at finite quark chemical potential. 
In Figl21 the dynamical quark mass is depicted as a function of the quark chemical 
potential. For = 313 MeV, the vacuum value for the dynamical quark mass 
is obtained. From Eq. (|2-9p at zero temperature, we obtain the relation between 
the quark number density pq and the quark chemical potential fiq. In Figj3l the 
quark number density is depicted as a function of the quark chemical potential. 
The vertical axis represents the quark number density divided by the normal nuclear 
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quark number density devided by normal nuclear density 




2 



340 



360 



380 



400 



quark chemical potential [MeV] 



Fig. 3. The quark number density versus quark chemical potential is depicted at T = 0. The 
dashed curve is not realized physically. 



Fig. 4. The pressure versus chemical potential is depicted at T = 0. The chiral phase transition 
occurs at fiq ~ 326 MeV. 

density, namely, Pq/p%- The dashed curves represent the unphysical region. Namely, 
the gap equation has multiple solutions in a certain region as is seen from Figl21 
Then, we must determine which solution is realized physically. For this purpose, the 
pressure is calculated from the thermodynamical potential density to in Eq. ()2-12p by 
the same subtraction method as the energy density per single nucleon in Eq. (|3-5p : 



where i = q for the quark matter. The physical solution corresponds to the one which 
gives the largest pressure. In FigU the pressure pq is depicted as a function of the 
quark chemical potential fig. From pg = 313 MeV to about 326 MeV, the low density 
solution is realized seen from Figj3j However, above pg w 326 MeV, the massless 
solution is physically realized. From FigiU it is seen that the chiral phase transition 
occurs at pg ^ 326 MeV. In this case, in the region from pg ~ 0.28/>^ to pg ~ 5Alp^ 
(from pb ~ 0.09/9^ to pB ~ l.SOp^ where pB represents the baryon number density) 
as is seen in FigjHl the first order chiral phase transition is realized and quark phases 



quark pressure [MeV/fm" ] 




-MT,p,) = mf^)){T,f,,) - mf'^W = 0,/i, = m,(r = 0))) 
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quark mass 
mq [MeV] 

~3DD 
200 
100 



320 325 

tlq [MeV] 

quark chemical potential 

Fig. 5. The solutions of the gap equations for quark mass are depicted as a function of the quark 
chemical potential /iq at T = with Gi^Ag = —68.4. 

Pq [/fm'^] 
2 

1.5 
1 

0.5 



340 360 380 400 420 
Hq [MeV] 

Fig. 6. The quark number density pq versus quark chemical potential Hq is depicted at T = with 
Gf„yl* — —68.4. The dashed curve is not realized physically. 

coexistence occurs. It seems that this density is rather small compared with the 
quark-hadron phase transition as is mentioned below. 

4.2.2. GivA^g = -68.4 

In the original NJL model with Giv = 0, it seems that the chiral phase tran- 
sition occurs at rather small density. One of possible reasons why the chiral phase 
transition occurs at rather small density may be that the strength of the attractive 
interaction between quark and antiquark is weak. So, the quark condensate is melt 
at rather small density. Thus, there exists one possibility that the attractive inter- 
action becomes strong by introducing the scalar-vector attractive interaction Giv for 
the original quark NJL model. 

However, there is no criterion to determine the value of Giv in this stage. Thus, 
as is the first attempt, we regard Gsv as a free parameter. First, we give the value 
G%A^g = -68.4, in which mg{pg/3 = p%) = 0.625mg{pg = 0). 

In FigEl the dynamical quark mass is depicted as a function of the quark chem- 
ical potential. For p,q = 313 MeV, the vacuum value for the dynamical quark mass 
is also obtained. In FiglHl the quark number density pg itself is depicted as a func- 
tion of the quark chemical potential p,q. The dashed curves represent the unphysical 
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Fig. 7. The pressure versus quark chemical potential is depicted and the chiral phase transition is 
shown at r = with Gj^A^ = -68.4. 

region, as is similar to FigjSl Namely, the gap equation has multiple solutions in 
a certain region. Then, we must determine which solution is realized physically by 
comparing with each pressure. In FigJTl the pressure for each branch is depicted as 
a function of the quark chemical potential /ig. 

From FigJTl it is seen that the chiral phase transition occurs at fiq ~ 326 MeV. 
In this case, in the region from pg ~ 0.38po to pg ~ 5.41/Jo (from ps ~ 0.12/)o to 
Pb ~ l.SOpo)) the chiral phase transition is realized and quark phases coexistence 
occurs. 

4.2.3. GivA^g = -225 

Finally, we give the value GivAg = —225 which leads to the dynamical quark 
mass at normal nuclear density as mg{pq/3 = p^) = 0.681mq(0). 

In FigiHl the dynamical quark mass is depicted as a function of the quark chem- 
ical potential. In FigjOl the quark number density pg is depicted as a function of the 
quark chemical potential pg. In FigfTOl the pressure of each branch is depicted as 
a function of the quark chemical potential. From Fig|10| it is seen that the phase 
with massless quark is realized in the region of pg > 392.5 MeV. In this case, at 

mq [MeV] 




340 360 380 

tlq [MeV] 

Fig. 8. The solutions of the gap equations for quark mass are depicted as a function of the quark 
chemical potential fiq at T = with Gf^/l** = —225. 
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|0q [/fm ] 



2 
1.5 

1 
0.5/ 




340 360 380 400 420 



Uq [MeV] 



Fig. 9. The quark number density pq versus quark chemical potential Hq is depicted at T = with 



Fig. 10. The pressure versus quark chemical potential is depicted and the chiral phase transition 
is shown at T = with GI„yl* = -225. 

Pq ~ 9.43/9^ {pB ~ 3.14/9^), the chiral phase transition may be realized. 



The main purpose of this paper is to investigate the quark-hadron phase tran- 
sition in the extended NJL model developed in this paper. In this section, we in- 
vestigate the realized phase by comparing the pressures of the nuclear and quark 
matters at zero temperature and the finite baryon chemical potential. It is shown 
that the first order quark-hadron phase transition is realized at zero temperature and 
the finite baryon chemical potential. In order to investigate the phase transition, the 
condition of the chemical equilibrium is demanded for the chemical potential pjy and 
Pq. Here, the chemical potential per one baryon at the same baryon density should 
be considered for the chemical equilibrium as follows: 



= -225. 




massless branch 



5. Quark-hadron phase transition 



Pn{T) = 3pq{T) . 




Also, the same pressures between the nuclear matter and the quark matters pq 
are necessary as 

Pn{T,pn) = Pq(T,3pq) . (5-2) 



14 



Y. Tsue, J. da Providencia, C. Providencia and M. Yamamura 



p.[MeV/fin] 




Sflq [MeV] 
jMn [MeV] 

Fig. 11. The quark-hadron phase transition is shown at T = with G^^A^ = —68.4. 

In (Pi/f(,) 




15 2 2.5 3 35 4 4.5 5 5.5 6 

biiryon number density / normal nuclear density /Jg / 

Fig. 12. The quark-hadron phase transition is shown at T = with Gj^A^ — —68.4. The vertical 
axis and the horizontal axis represent the pressure \n{pi/po) and the baryon number density 
divided by the normal nuclear density Pb/p%, respectively. Here, we use po = 1.0 MeV/fm^. 

The pressure has aheady been defined in Eq. ()4-2p . 

The crossing point of each pressure for nuclear matter and the quark matter 
presents the coexistence point of the nuclear and the quark phases and the phase 
with larger value of the pressure is realized except for the crossing point. 

In FigHTJ the pressures for the nuclear matter and the quark matter are depicted 
as functions of the nuclear chemical potential /Utv and the quark chemical potential 
multiplied by 3, 3fiq, respectively, in the case Giv^q = —68.4. In this figure, it is 
shown that the quark-hadron phase transition occurs at a critical chemical potential 
value fic which is about Hci= A* at = w 1236 MeV. From this figure, the nuclear or 
hadron phase is realized in the small chemical potential region fi]\f < fic- However, 
in the large chemical potential region 3fj,q > fic, the quark phase is realized. In 
terms of the baryon number density by using the relation of the density and 
chemical potential shown in FigEl from pq ~ 10.9po {PB = Pq/3 ~ 3.Q3po) to 
Pn{= Pb) ~ 2.65p0) the nuclear and quark phases coexist and the first order quark- 
hadron phase transition is realized. This situation is depicted in Fig|12[ 

For the other parameter values for Gf^, the same figure is obtained. The reason 
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is as follows: For these model parameters used in this paper, that is, GivA^ = 0, 
—68.4 and —225, the quark-hadron phase transition point is in the chiral symmetric 
phase, that is, niq = and {ipqipq) = 0. Namely, the quark-hadron phase transition 
from hadron phase to quark phase occurs after the chiral phase transition from chiral 
broken phase to chiral symmetric phase in these model parameters. This behavior 
for the phase transitions is also seen in Ref{T2j). Then, the pressure of quark matter 
Pq in Eq. ()4-2p with T = does not depend on Giv because (T-Lq^^) in Eq. (|2-13p and 
fj,q in Eq. (j2-6|) does not depend on Gs„ due to {i^qipq) = 0. Thus, the behavior seen 
in FigfT2] is not changed. 

§6. Summary and concluding remarks 

In this paper, the quark-hadron phase transition at finite baryon chemical po- 
tential has been described in the extended NJL model, in which both the nucleon in 
the nuclear matter and the quark in the quark matter were treated as fundamental 
fermions with a number of color, Nc, being one for nucleon and being three for quark. 
In this paper, as the first attempt of the investigation for the quark-hadron phase 
transition in the extended NJL model, we only dealt with the symmetric nuclear 
matter and the quark matter without quark-quark correlation leading to the pairing 
instability. 

As for the symmetric nuclear matter, the saturation property has been well 
reproduced in this model. As for the quark matter, there is one free parameter, 
that is, the coupling strength of the scalar- vector eight-point interaction, This 
model parameter controlls the chiral phase transition point and/or the strength of 
the partial restoration of the chiral symmetry in the nuclear medium. In this paper, 
we did not fix its value by using other physical quantity. This is a future problem. 

We have calculated the pressures of nuclear matter and quark matter, respec- 
tively. Then, by comparing the pressure of nuclear matter with that of quark matter, 
the realized phase is determined. As a result, the first order quark-hadron phase tran- 
sition is obtained at finite density and zero temperature. It should be noted here 
that, for the adopted parameter Giv used in this paper, the deconfinement phase 
transition occurs after the chiral symmetry restoration in the nuclear matter. Thus, 
the phase transitions in the quark phases may be hidden because they occur before 
quark deconfinement. 

It is interesting to investigate the quark-hadron phase transition at finite temper- 
ature and baryon chemical potential because the phase diagram in the QCD world 
is not understood completely at present. It is also interesting to study the phase 
transition between the neutron matter and quark matter at finite density, which 
has important implications to the physics of neutron star. Of course, on the side of 
quark phase, the color superconducting phase should be taken into account. These 
are future problems in this model calculation. 
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